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ABSTRACT 

We present angular power spectra of the cosmic microwave background radiation anisotropy 
due to fluctuations of the Sunyaev-Zel'dovich (SZ) effect through clusters of galaxies. A 
contribution from the correlation among clusters is especially focused on, which has been 
neglected in the previous analyses. Employing the evolving linear bias factor based on the 
Press-Schechter formalism, we find that the clustering contribution amounts to 20 — 30% of the 
Poissonian one at degree angular scales. If we exclude clusters in the local universe, it even 
exceeds the Poissonian noise, and makes dominant contribution to the angular power spectrum. 
As a concrete example, we demonstrate the subtraction of the ROSAT X-ray flux- limited 
cluster samples. It indicates that we should include the clustering effect in the analysis of the 
SZ fluctuations. We further find that the degree scale spectra essentially depend upon the 
normalization of the density fluctuations, i.e., as, and the gas mass fraction of the cluster, rather 
than the density parameter of the universe and details of cluster evolution models. Our results 
show that the SZ fluctuations at the degree scale will provide a possible measure of erg, while 
the arc-minute spectra a probe of the cluster evolution. In addition, the clustering spectrum 
will give us valuable information on the bias at high redshift, if we can detect it by removing 
X-ray luminous clusters. 



Subject headings: cosmology: cosmic microwave background - galaxies: clusters 



1. Introduction 

The cosmic microwave background radiation (CMBR) anisotropy is now known as the most powerful 
probe of our Universe (e.g. Hu, Sugiyama & Silk 1997). In addition to the primary anisotropy which will give 
us clues to most of the cosmological parameters within 10% accuracy (Bond, Efstathiou & Tegmark 1997), 
there is the secondary anisotropy caused by several processes in the low-redshift universe, e.g., non- linear 
evolution of gravitational structures (Rees & Sciama 1968), peculiar velocity field of re- ionized clumps 
(Ostriker & Vishniac 1987), and the Sunyaev-Zel'dovich (SZ) effect through clusters of galaxies. The SZ 
effect is the brightness change of the CMBR due to inverse Compton scattering by hot electrons (Zel'dovich 
& Sunyaev 1969; Sunyaev & Zel'dovich 1972). So far the brightness decrement in the Rayleigh- Jeans band 
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has been detected toward tens of clusters (see, Birkinshaw 1999 for review), and we have recently detected 
the increment in the Wien band (Komatsu et al. 1999). Thanks to such specific spectral features, one will 
be able to separate the SZ fluctuations (SZF) rather accurately from any other signals using future satellite 
observations, e.g., that of PLANCK sm^reyoi (Hobson et al. 1998). Furthermore, the SZF should provide 
a probe of high-z clusters yet to be detected otherwise, because the SZ brightness is free from the redshift 
diminishing effect (Birkinshaw 1999) unlike the X-ray brightness. 

Although several theoretical studies have been done on the SZF (Rephaeli 1981; Ostriker & Vishniac 
1986; Cole & Kaiser 1988; Schaeffer & Silk 1988; Markevitch et al. 1991; Makino & Suto 1993; Bartlett 
& Silk 1994; Persi et al. 1995; Bond & Myers 1996; Atrio-Barandera & Miicket 1999), they mainly focus 
on the Poissonian component only and much less attention is paid to that arisen from the cluster-cluster 
correlation. In this Letter, we point out that the cluster-cluster correlation can in fact make a non-negligible 
contribution to the SZF at degree scales and explicitly predict its angular power spectra. We further 
investigate possibilities of exploring the nature of density fluctuations, cluster gas distributions, and 
clustering properties of high-redshift clusters, by means of the SZF. 



2. Model 
2.1. Angular Profile of a Cluster 

The gas density profile of each cluster is approximated by the spherical isothermal /3-model, 

-3/3/2 

, where is the core radius, and Pgaso is the central gas mass density. 



1 + {r/r,f 



Pgas[r) = PgasO 

To account for the fact that a real cluster has a finite size, we introduce the Gaussian filter on the edge 
at radius R, i.e., Pgasir) Pga,s{r)e^^ , where ^ is a fudge factor to normalize the total gas mass of a 
cluster. Assuming /3 = 2/3 for simplicity, the gas mass contained within R is 

Mgas = / dVpg_^s{r) = 47rpgasore {p - arctanp) , (1) 
Jo 

where p = R/r^. On the other hand, 

o 

dFpgas(r)e-'-'/«^' = Anp^^orlHp), (2) 



where 

h{p) ^ p^i - |ei/«^'Erfc (VW) , (3) 

and Erfc(?i) is the complementary error function. Matching equations (|l|) and (^ at p 3> 1 gives — A/it. 

The above density profile yields the angular temperature profile projected on the sky due to the SZ 
effect in terms of the "Compton j/-parameter" (Zel'dovich & Sunyaev 1969) as 

AT(6) 



^cmbr 



9{x)y{e), (4) 




y{e) = arneorc f^^) x =Erfc ( ^/ ^ ' ^"J.^""' ) , (5) 

where Tcmbr — 2.728K (Fixsen et al. 1996), g{x) = a;coth(x/2) — 4, a; = hi/ / kBTcmhr , n-eo is the central 
electron number density, Te is the electron temperature, and other physical constants have their usual 
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meanings. The angular core size 9c is related to Tc via 9c — rc/dA{z), where (Ia is the angular diameter 
distance to a cluster at redshift z. 

We model that clusters have virialized at z and relate the quantities R, Tg, rigo and Tc to the total 
mass M and redshift z based on the spherical infall model (Peebles 1980, §19). We assume that R 
is equal to the virial radius, R — [3M/47ryOuniv(2)Ac(z)]^/'^, where Puniv{z) is the mean density of the 
universe, and Ac{z) is the mean density of a cluster in units of Puniv(-z) given in Lacey & Cole (1993) 
and Nakamura & Suto (1997) for an arbitrary cosmology. Then the electron temperature is expressed as 
fceTe = rj^nipGM /iR, using the mean molecular weight p, and a fudge factor of order unity ry, for which 
we adopt fi — 0.59 and rj = j3~^ = 3/2, respectively. The central electron number density is derived with 
the help of equation (||) as rieo = {Pf,&sol'mp){X + \)/2— (A//gas/?Tip47rrj?/i(p)) {X + l)/2, where we fix the 
hydrogen abundance a± X = 0.76 and assume that the gas mass fraction is equal to the cosmological mean, 
/gas = Afgas/Af = r^b/^^o (White et al. 1993). Given large uncertainties in the formation of cluster core 
regions, we analyze two different scenarios for the evolution of r^, i.e., the self- similar (SS) model (Kaiser 
1986), and the entropy-driven (ED) model (Kaiser 1991; Evrard & Henry 1991; Bower 1997). In the SS 
model, Tc is simply proportional to the virial radius, 

Pss{M,z) = constant. (6) 

In the ED model, on the other hand, the core is in a minimum entropy phase, Smin = cy ln(re/Pgaso)j ''^^'^ 
its evolution is specified by Sm\n(z) = Smin(-z = 0) + cyeln(l + z) with a free parameter e, yielding 

1/2 



Pcd(M,z) 



3pgaso(M, z)//ga 
Puniv(^:)Ac(z) 

TI^{M,z) 



AV {z) 

where we have used M — (47r/3)puniv(^) Ac(z)i?'^ ~ 47r(/9gaso//gas)?'cP for p 1. The adiabatic index is 
fixed at 7 = 5/3 throughout. Since Mushotzky & Scharf (1997) obtained the constraint of e = ± 0.9 from 
X-ray cluster samples, we demonstrate three cases of e = —1, 0, 1 in the same way as Atrio-Barandela & 
Miicket (1999). We normalize the proportional constants in both models at M = 10^^h~^MQ and z — Q 
assuming r c{lO^^ h'^ Mq, z ^ 0) = 0.15/i~iMpc. 



2.2. Angular Power Spectrum 

One conventionally expands the angular two-point correlation function of the temperature distribution 
in the sky into the Legendre Polynominal: 

/ AT AT \ 1 

/ -(n)- {n + e)) = —yCiPi{cos9). (8) 



\ -^cmbr ^cmbr / ^TT ^ 

iP) IC) (P) 

Since we are considering discrete sources, we can write Ci ~ + Cl , where C; is the contribution 
from the Poissonian noise and C\ is from the correlation among clusters (Peebles 1980, §41). For 
convenience, we define the power spectrum of j/-parameter, C^^ = Cijg^ix)^ which is independent of 
frequency. Following Cole & Kaiser (1988) and employing Limber's equation in the cosmological context 
(e.g., Peebles 1980, §56), we get an integral expression of (72/2/(^,0). 



dz ./m '^-^ 
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where V{z) and r{z) are the comoving volume and the comoving distance, respectively, <^i{M,z) is the 
"biased temperature function" : 



<^i{M,z) 



_ dn(M, z) 
" dM 



h{M,z)yi{M,z), 
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and yi is the angular Fourier transform of y{9). In equations (^) and (|lO|), Zdoc is the redshift of photon 
decoupling, and we choose M^in = 5 x lQ^^h~^MQ and M^ax = 5 x IQ^^H'^Mq to bracket the mass scale 
of clusters. The mass function dn/dM is computed using the Press-Schechter formula (Press & Schechter 
1974): 



dn{M, z) 
dM 
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where a{M) is the present-day rms mass fluctuation within the top-hat filter, v{M,z) = 5c{z) / D{z)a{M) 
is the "peak-height threshold", D{z) is the linear growth factor of density fluctuations, and 5c{z) is the 
threshold overdensity of spherical collapse given in Lacey & Cole (1993) and Nakamura & Suto (1997) 
for an arbitrary cosmological model. We assume that the matter power spectrum, P-ai(k,z), is related to 
the power spectrum of the cluster correlation function, Pc{k, Mi, M2, z), via the time-dependent linear 
bias factor, b{M,z), i.e., Pc{k, Mi, M2, z) = h{Mi, z)h{M2, z)D'^{z)P,a{k, z = 0). We adopt for b{M,z) an 
analytic formula derived by Mo & White (1996): b{M,z) = 1-1- (^iy^{M,z) — l) /Sc{z), which shows an 
excellent agreement with numerical simulations at the cluster mass scale (Jing 1998, 1999). To calculate 
Pm(fc), we employ the Harrison-Zel'dovich spectrum and the transfer function of Bardeen et al. (1986) with 
the shape parameter given in Sugiyama (1995). Unless stated explicitly, we assume the ACDM cosmology 
with JIq = 0.3, Ao — 0.7, fib — 0.05, h = 0.7 and as = 1.0, which agrees with both the COBE normalization 
and the cluster abundance (e.g. Kitayama & Suto 1997). 



3. Results and Discussions 

3.1. Poissonian and Clustering Contributions to the Sunyaev-Zel'dovich Fluctuations 

Figure |^ shows the predicted power spectra of y-parameter, Cf^^-^'^\ in the SS and three ED 
(e = —1,0,1) models. The correlation term Cf^^'^"^ amounts to around one-fifth of the Poissonian term 
(jvv(p) ^-(- I < j^gg^ gj^^ show later that Cf^^^^ even dominates Cf^'^-^^ in some practical situations 

(Figure I). 

For comparison, the primary anisotropy expected in units of the Rayleigh- Jeans temperature is overlaid 
on the same figure. The SZF signals are shown to be sub-dominant except at I ^ 2000. However, the 
simulations of Hobson et al. (1998) indicate that the SZF can be separately reconstructed from the CMBR 
map by the PL^A'^C/ir satellite with great accuracy especially at large angular scales (/ ^ 500). This is fully 
due to the specific spectral features of the SZ effect in contrast to the primary signal and other secondary 
sources. The angular resolution of the PLANCK, on the other hand, seriously limits the accuracy of 
reconstruction at small angular scales (/ ^ 1000). 

The SZF spectra at I ^ 500 is quite insensitive to the details of core evolution models. This feature can 
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be understood by looking at the angular Fourier transformation of y{0) (see eq. ||^): 



n^QT^e^e"^'^-' /l {Wc > 1), 



The spectrum at large angular scales is solely determined by the gas mass: Ui^g-^ oa A/gas, while the small 
scale spectrum depends on the core radius: Ui^g-i oc Mg^^e~^^'' / {Itc). In what follows, therefore, we focus 
mainly on the large scale spectrum and consider only the SS model. 

In practical observations, one can remove from the reconstructed SZF map bright clusters already 
detected in X-ray observations, e.g., by iJOS'^T satellite. Figure || illustrates the expected SZF spectra after 
excluding the clusters above some threshold X-ray fluxes, S'jf (0.5 — 2keV) — 10~^^ and 10""'^'^ erg cm~^ s~^. 
These values roughly correspond to the limiting fluxes in the iZ 05*^4 T" All-Sky Survey (Ebeling et al. 1997) 
and the ROSAT De&p Cluster Survey (Rosati et al. 1998). To compute the X-ray fluxes of clusters, we 
have used the same method as described in Kitayama & Suto (1997), assuming the following form for the 
Lx ~ T relation: 

_fcBT V 

with igkcv = 2.9 X lO^'^/i-^, a = 3.4 and C = (David et al. 1993; Mushotzky & Scharf 1997), and the 
X-ray band correction according to the emission model by Masai (1984). 



= iekcv ( 7r|-^ j (1 + 2)"^ ergs-i, (14) 



We find that the Poissonian term C;^^*^'^-' is greatly reduced by the subtraction of X- 
ray bright clusters, and the correlation term C\^'^'~^^ indeed dominates at I < 200. This 
reflects the greater contribution of the correlation term Cf^^^^ at higher z and smaller I, since 
Cfy(c)/Cfy(P) - n{M,z)Pc{k = l/r{z)) - n{M, z)P{M, z)D^{z)r{z)/l, where P,„(fc) - fc-^ on cluster 
scales. Therefore, we should take into account C^^'^'-^^ in the analysis of the SZF, and it will possibly be 
detected by the PLANCK ^Nith removing the 05^ T clusters. 

Now, what does Cf^^Q^Q mean after this subtraction ? It should represent the clustering properties of 
high-z clusters currently unresolved in X-ray observations. It should be noted that the above predictions 
depend on the linear bias assumption and the evolution history of the bias factor. Thus, if we could measure 
(jvv{C) separately from Cf^^'P\ it would in turn provide the information on the bias at high- 2;. 



3.2. Implications on Cosmology 

Figure || further illustrates the as dependence of C^jj^p'^'' in three representative cosmological models: 
ACDM, SCDM (f7o 1, A 0, = 0.5) and OCDM {Uo = 0.45, A 0, = 0.7). The SZF spectra are 

(C) 

extremely sensitive to erg while less sensitive to Qq and Aq. The correlation term is more sensitive to 
(78 than Cl^^ , and can be fitted roug hly by C/^) oc al'^cl^^ Since the PLANCK wiW provide the accurate 
angular power spectrum of the SZF at the degree angular scales (Hobson et al. 1998), they will provide a 
plausible measure of erg. Furthermore, figures |l| and ^ suggest that the degree scale spectrum is free from 
both unknown core physics of clusters and density parameters of the universe. 

Bond, Efstathiou & Tegmark (1997) found that the expected la error of as determined by the 
PLANCK using the primary anisotropies alone is as large as 20%, which corresponds to changes in C; by 
a factor of 3 and C^'^'^ by a factor of 4. Figure ^ shows the total (Poissonian -I- clustering) and Poissonian 
angular power spectra of the SZF in units of the Rayleigh- Jeans temperature, ATi = Tcmbr \/^(H-l)C/ . We 
show both cases of including and excluding the X-ray bright clusters with Sx > 10"^'^ erg cm^^ s^^. It 
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illustrates how much the SZF change within 20% error of erg. Notice that C; without X-ray clusters has 

(c) (P) 

a bending structure at I ^ 200, below which C, exceeds . The spectral shape also changes with 
as because and depend differently on cjg. Such features should be quite helpful in detecting 

the clustering contribution separately from the Poissonian one and probing the structure formation in the 
high-2 universe. For comparison, we also plot the primary anisotropy spectrum and observational upper 
bounds obtained by OVRO (Readhead et al. 1989) and ATCA (Subrahmanyan et al. 1998). 

The SZF should play a complementary role to the primary anisotropics in the determination of as- 
and the clustering contribution hopefully brings us fruitful information of high-z bias. It should be noticed, 
however, that these results are based upon the assumption of /gas = Ob/fio- We found a rough scaling 
relation: C^^^ oc f'^^^a^Vt^h, at degree angular scale, and it gives C\^^ oc fig when we specify /gas without 
cosmological parameters. Uncertainties in /gas yield systematic errors in the amplitude of predicted SZF 
angular spectrum. These points will be discussed in greater detail in future publications. 

We thank Naoshi Sugiyama and Yasushi Suto for valuable suggestions, Toshifumi Futamase, Makoto 
Hattori, Masahiro Takada and Keiichi Umetsu for frequent discussions and critical comments. We 
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REFERENCES 

Atrio-Barandela, F., & Miicket, J. P. 1999, ApJ, 515, 465 

Bardeen, J. M., Bond, J. R., Kaiser, N., & Szalay, A. S. 1986, ApJ, 304, 15 

Bartlett, J. G., & Silk, J. 1994, ApJ, 423, 12 

Birkinshaw, M. 1999, Phys. Rep., 310, 97 

Bond, J. R., Efstathiou, G., & Tegmark, M. 1997, MNRAS, 291, L33 
Bond, J. R., & Myers, S. T. 1996, ApJS, 103, 1 
Bower, R. G. 1997, MNRAS, 288, 355 

Catelan, P., Lucchin, F., Matarrese, S., & Porciani, C. 1998, MNRAS, 297, 692 

Cole, S., & Kaiser, N. 1988, MNRAS, 233, 637 

David, L. P., Slyz, A., Jones, C, Forman, W., Vrtilek, S. D., & Arnaud, K. A. 1993, ApJ, 412, 479 
Ebcling, H., Edge, A. C, Fabian, A. C, Allen, S. W., & Crawford C. S. 1997, ApJ, 479, LlOl 
Evrard, A. E., & Henry, J. P. 1991, ApJ, 383, 95 

Fixsen, D. J., Cheng, E. S., Gales, J. M. Mather, J. C, Shafer, R. A., & Wright, E. L. 1996, ApJ, 

473, 576 

Hobson, M. P., Jones, A. W., Lasenby, A. N., & Bouchet, F. R. 1998, MNRAS, 300, 1 
Hu, W., Sugiyama, N., & Silk, J. 1997, Nature, 386, 37 
Jing, Y. P. 1998, ApJ, 503, L9 
Jing, Y. P. 1999, ApJ, 515, LI 



- 7- 



Kaiser, N. 1986, MNRAS, 222, 323 

Kaiser, N. 1991, ApJ, 383, 104 

Kitayama, T., & Suto, Y. 1997, ApJ, 490, 557 

Komatsu, E., Kitayama, T., Suto, Y., Hattori, M., Kawabe, R., Matsuo, H., Schindler, S., & 
Yoshikawa, K. 1999, ApJ, 516, LI 

Lacey, C, & Cole, S. 1993, MNRAS, 262, 627 

Makino, N., & Suto, Y. 1993, ApJ, 405, 1 

Markevitch, M., Blumenthal, G. R., Forman, W., Jones, C, & Sunyaev, R. A. 1991, ApJ, 378, L33 

Masai, K. 1984, Ap&SS, 98, 367 

Mo, H. J., & White, D. M. 1996, MNRAS, 282, 347 

Mushotzky, R. F., & Scharf, C. A. 1997, ApJ, 482, L13 

Nakamura, T. T., & Suto, Y. 1997, Prog. Theor. Phys., 97, 49 

Ostriker, J. P., & Vishniac, E. T. 1986, ApJ, 306, L51 

Peebles, P. J. E. 1980, The Large Scale Structure of the Universe (Princeton: Princeton Univ. Press) 
Persi, F. M., Spergel, D. N., Cen, R., & Ostriker, J. P. 1995, ApJ, 442, 1 
Press, W. H., & Schechter, P. 1974, ApJ, 187, 425 

Readhcad, A. C. S., Lawrence, C. R., Myers, S. T., Sargent, W. L. W., Hardebeck, H. E., & Moffet, 
A. T. 1989, ApJ, 346, 566 

Rees, M. J., & Sciama, D. W. 1968, Nature, 517, 611 

Rosati, P., Delia Ceca, R., Norman, C, & Giacconi, R. 1998, ApJ, 492, L21 
Shaeffer, R., & Silk, J. 1988, ApJ, 333, 509 

Subrahmanyan, R., Kesteven, M. J., Ekers, R. D., Sinclair, M., & Silk, J. 1998, MNRAS, 298, 1189 
Sugiyama, N. 1995, ApJS, 100, 281 

Sunyaev, R. A., & Zel'dovich, Ya. B. 1972, Comments Astrophys. Space Phys., 4, 173 
White, S. D. M., Navarro, J. F., Evrard, A. E., & Frenk, C. S. 1993, Nature, 366, 429 
Zel'dovich, Ya. B., & Sunyaev, R. A. 1969, Ap&SS, 4, 301 



- 8 - 




Fig. 1. — Angular power spectra of the Sunyaev-Zel'dovich fluctuations in terms of y-parameter (see text 
for definitions). The Poissonian and clustering spectra are shown for different core evolution models. Solid 
lines represent the self-similar model, while others the entropy-driven model with e = —1,0, 1 from top to 
bottom. Also plotted for reference is the primary temperature anisotropy expected in the Rayleigh-Jeans 
band {right vertical axis). 
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Fig. 2. — Effects of subtracting the clusters already resolved in the X-ray band. Filled circles show the 
angular power spectra without the subtraction. Dashed and solid lines are the spectra after removing the 
X-ray clusters above the flux limit of 10~^^ and 10^^'^ erg cm~^ s""'^ in the 0.5 — 2keV band, respectively. 
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Fig. 3. — Dependences on erg of angular power spectra at / = 100 in three representative cosmological models, 
ACDM (solid lines), OCDM (dashed) and SCDM (dotted). The upper panel shows the Poissonian (thin lines) 
and clustering (thick lines) power spectra separately, while the lower panel illustrates ratios between them. 
Filled triangles indicate reference lines oc Cg and CTg"^'^ in the upper and lower panels, respectively. Filled 
circles in the lower panel represent the tig value inferred from the cluster abundance in each cosmological 
model (Kitayama & Suto 1997). 
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Fig. 4. — Total (thick lines) and Poissonian (thin lines) angular power spectra of the Sunyaev-Zel'dovich 
fluctuations in the Rayleigh-Jeans band excluding the X-ray bright clusters with Sx{0-5 — 2keV) > 
10^^"^ erg cm"^ s~^. Dotted, solid and dashed lines represent the cases of erg — 0.8, 1.0 and 1.2, respectively. 
Circles show the total power spectrum with the X-ray bright clusters for erg = 1.2. Also plotted for reference 
are the expected primary anisotropy and the observed upper limits reported by OVRO (20 GHz) and ATCA 
(8.7 GHz). 



